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The surface tension of cold and dense QCD phase transitions has appeared recently as a key 
ingredient in different astrophysical scenarios, ranging from core-colapse supernovae explosions 
to compact star structure. If the surface tension is low enough, observable consequences are 
possible. Its value is however not known from first-principle methods in QCD, calling for effective 
approaches. Working within the framework of homogeneous nucleation by Langer, we discuss 
the steps that are needed to obtain the nucleation parameters from a given effective potential. As a 
model for deriving the effective potential for the chiral transition, we adopt the linear sigma model 
with constituent quarks at very low temperatures, which provides an effective description for the 
thermodynamics of the strong interaction in cold and dense matter, and predict a surface tension 
of £ ~ 5-15 MeV/fm 2 , well below previous estimates. Including temperature effects and vacuum 
logarithmic corrections, we find a clear competition between these features in characterizing the 
dynamics of the chiral phase conversion. 
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One of the many faces of QCD that has been attracting continuous attention in the last years is 
how the matter dictated by the fundamental theory of strong interactions behaves at extreme con- 
ditions of temperature and/or density. The fact that the energy scales for the QCD phase transitions 
are currently under reach of different heavy-ion colliders creates an intense demand of theoretical 
and phenomenological development. Over the past decade RHIC has established the formation of 
a new state of matter, characterized by a clear fluid behavior at the partonic level, differing from 
the predicted picture of the naive weakly-interacting quark gluon plasma [jlj]. LHC stretches the 
energy frontier also for heavy-ion collisions, bringing up the statistics and allowing for a careful 
study of fluctuations that might reveal the features of the initial state of the collisions. 

The richest structure in the high energy QCD phase diagram seems however to be deep within 
the region with an asymmetry between the particle and antiparticle content of the medium. In this 
direction many experimental efforts are dedicated to create matter with a finite baryon-chemical 
potential: the aim is the search for the phase boundaries in the QCD phase diagram and its proper- 
ties. The accelerators however cannot reach the domain of very low temperature and high chemical 
potential, which is also the region where theorists have less guidance for what QCD predicts, since 
first-principle lattice calculations are not possible [J2|]. 

Astrophysical phenomena might however shed some light in this partially obscure region of 
the QCD phase diagram. Extremely high densities and relatively low temperatures are believed 
to exist in the core of ultracompact objects, so that it is probable that the QCD phase transitions 
play an important role in the structure and dynamics of such systems. Even though in general 
observations are much less controlable than laboratory experiments, the astronomical measurement 
techniques and instrumentation have much developed in the last decade and this might be the 
dawn of a precision era in which astrophysics really constrains nonperturbative QCD phenomena. 
One concrete example of astronomical observation that has reached the capability of constraining 
models for strong interactions in the nonperturbative regime is that of Ref. ^ in which the error 
bars of less than 5% in the 2 solar masses measurement rule out different equations of state. 

The surface tension between partonic and hadronic phases of QCD seems to be a key quan- 
tity in connecting microscopic modeling of cold and dense QCD with astrophysical applications: 
recently it has been put forward that a low enough QCD surface tension should allow for possi- 
ble observable consequences in different phenomena related to compact objects. It was shown for 
instance that deconfinement can happen during the early post-bounce accretion stage of a core col- 
lapse supernova event, which could result not only in a delayed explosion but also in a neutrino 
burst that could provide a signal of the presence of quark matter in compact stars M\. However, 
as was discussed in detail in Ref. [||] (see also Ref. [||] ) those possibilities depend on the actual 
dynamics of phase conversion, more specifically on the time scales that emerge. In a first-order 
phase transition, as is expected to be the case in QCD at very low temperatures, the process is 
guided by bubble nucleation (usually slow) or spinodal decomposition ("explosive" due to the van- 
ishing barrier), depending on how fast the system reaches the spinodal instability as compared to 
the nucleation rate []7||. Nucleation in relatively high-density, cold strongly interacting matter, with 
chemical potential of the order of the temperature, can also play an important role in the scenario 
proposed in Ref. [j|], where a second (little) inflation at the time of the primordial quark-hadron 
transition could account for the dilution of an initially high ratio of baryon to photon numbers. 
Moreover, significantly different compact star structures are obtained if one considers the possibil- 
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ity of a layer of a quark-hadron mixed phase in the core of compact stars [Qj. A key ingredient in 
all these scenarios is, of course, the surface tension, since it represents the price in energy one has 
to pay for the mere existence of an interface between quark and hadron phases. If the cost in energy 
is too high, quark-hadron mixed phases will not be favourable and similarly nucleation time scales 
will be too long for these astrophysical phenomena to take place. 

Computing the surface tension for the QCD transitions with first-principle methods is a pro- 
hibitively complicated task. It is a genuine nonperturbative quantity, in the sense that it cannot be 
obtained by a naive expansion around one minimum of the effective action, and is also inaccessible 
by current first-principle lattice QCD simulations, due to the sign problem. One has therefore to 
resort to effective models. Estimates for the surface tension between a quark phase and hadron mat- 
ter were considered previously in different contexts. In a study of the minimal interface between 



a color-flavor locked phase and nuclear matter in a first order transition, the authors of Ref. [ |10| ] 
use dimensional analysis and obtain £ ~ 300 MeV/fm 2 assuming that the transition occurs within 
a fermi in thickness. Taking into account the effects from charge screening and structured mixed 



phases, the authors of Ref. [11] provide estimates in the range of 50 — 150 MeV/fm 2 but do not 
exclude smaller or larger values. These values of the surface tension are too high, rendering the 
nucleation of quark matter bubbles improbable in the context of compact objects. 



We have estimated [12] the surface tension for the cold and dense chiral phase transition within 



a well-established chiral model (namely, the linear sigma model with constituent quarks (LSMq)) 



and assuming the homogeneous nucleation scenario [13]. Our results indicate that the surface 
tension in the cold and dense regime can be considerably smaller than previous estimates, possibly 
allowing for the interesting observable astrophysical phenomena described above. 

This paper is organized as follows. In Section [IJ, we delineate how, within the framework 
of homogeneous nucleation, it is possible to derive the nucleation parameters starting from an 
effective potential for the chiral order parameter V e s(c). Section || contains a short presentation of 
the effective model we adopted for the QCD chiral phase transition and presents the results for the 
effective potential V e ff (a) for three different cases. Our results for the surface tension are discussed 
in Section || and Section |] contains final remarks. 

1. How to estimate the surface tension? 

As stated above, with the aim of quantifying the surface tension associated with the cold and 
dense chiral phase transition in QCD, we will adopt the framework of homogeneous nucleation 
proposed by Langer JOt). It is based on the construction of a coarse-grained free energy functional 
&[o\ for the configuration of the order parameter field a from a given effective potential V e ff as 
follows: 



&[o] = An fr 2 dr 



i(£i — 



(i.i) 



where spherical symmetry has been assumed. 

Around the domain of temperature T and chemical potential /I in which a first-order phase 
transition takes place, a metastable phase appears. More precisely, metastability manifests as the 
existence of (at least) two minima of the effective potential with a maximum in between inside the 
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T\- spinodal lines 
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Figure 1: Cartoon of the temperature versus chemical potential phase diagram with a first order critical line 
(solid) and the spinodal curves (dashed) defining the metastable region. 



region defined by the two spinodal lines, as depicted in Fig. [j]. Within the so-called metastable 
region, the coarse-grained free energy functional in Eq. ( |l.l| ) develops a nonuniform extremum 
configuration of the bubble type tpb(r;R = R c ), connecting the metastable minimum to the absolute 
one (for an illustration cf. Fig. |2|). This unstable nonuniform extremum is called the critical bubble 
and is intimately related to the dynamics of phase conversion via bubble nucleation. Differently 
from the uniform configuration sitting at the absolute minimum, the critical bubble is unstable, 
corresponding actually to a maximum of the free-energy & ' . 

If the system goes through an out-of-equilibrium stage in which the external parameters T and 
/I are modified fast enough as compared to the time scales of the microscopic dynamics, then the 
system may get trapped in the metastable phase. In this scenario, fluctuations (brought about by 
the interaction with a thermal or particle reservoir) may generate bubbles of the absolute minimum 
inside the metastable medium. These fluctuation bubbles will then evolve according to the com- 
petition between volume and surface terms, as dictated by the free-energy difference between the 
configurations with and without a bubble of radius R: 



Ayr 

A^ = ^[<j> b {r,R)]-^[<l> b (r,R = 0)} = -—R 3 AP + 47iR 2 T l , 



(1.2) 



whose general form is sketched in Fig. ||. The volume term encodes the gain in energy resulting 
from the conversion to the true equilibrium state and is proportional to the pressure difference 
AP > between the phases. On the other hand, the surface term accounts for the energy spent in 
the construction of a surface between the phases, being given essentially by the surface tension £. 
As illustrated in Fig. |3|, the free-energy difference AJ^~ shows a maximum associated with 
the critical bubble (R = R c ) and fluctuation-generated configurations containing bubbles with radii 
R < R c will tend to the uniform metastable vacuum in order to minimize the energy. For radii bigger 
than the critical radius, the minimization of energy implies that the droplet will grow and eventually 
complete the phase conversion. The critical radius goes to zero at the spinodal curves, where the 
metastable false vacuum becomes unstable (the barrier disappears) and the phase conversion occurs 
explosively via the spinodal decomposition process. 



Finally, from Eq (1.2) it is clear that the surface tension can be obtained via the knowledge 
of the effective potential V e fp and the critical bubble profile <pt{r',R = R c )- Since the critical bubble 
itself is an unstable extremum of the coarse-grained action J? in Eq. ( |1.1| ), it is ultimately obtained 
from the effective potential V e ff as well. In the thin wall approximation (valid when the critical 
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Figure 2: Example of bubble profile. 
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Figure 3: Sketch of the difference of coarse- 
grained free energy between the configura- 
tions with and without a bubble of radius R. 



bubble has a thin interface region as compared to its radius, i.e. % <C R c ', cf. Fig. ^) for a quartic 
effective potential, the solution <j)b{r',R = R c ) can be written analytically in terms of the Taylor 
coefficients of the potential. 



2. Effective model for the chiral phase transition 

In order to estimate the chiral surface tension for cold and dense QCD matter, it suffices there- 
fore, within the framework of homogeneous nucleation presented above, to provide an effective 
potential Veff(cr) describing the chiral phase transition for low temperatures and high chemical 
potentials. 

Here we consider the linear sigma model coupled with constituent quarks (LSMq) at zero or 
low temperature and finite quark chemical potential as a model for the thermodynamics of strong 
interactions in cold and dense matter. Since we are concerned with the chiral phase transition, we 



neglect the pions [ 14] and work with the following Lagrangian: 



1 



where 



if = y f [i'fd il -m q -go]yf + -d il cd ll o-U{o), 



U{c) = ^{a 2 -v 2 ) 2 -ha 



(2.1) 



(2.2) 



is the self-interaction potential for the sigma meson, exhibiting both spontaneous and explicit 
breaking of chiral symmetry. The massive fermion fields \\ff (f = 1, ■ •• ,Nf = 2) represent the 
up and down constituent-quark fields. The scalar field a plays the role of an approximate order 
parameter for the chiral transition, being an exact order parameter for massless quarks and pions. 
We compute the effective potential 



v efl ((j) = u(cj)+a™((j) 



(2.3) 



for the sigma condensate integrating over the quark fields and keeping quadratic fluctuations of 
the sigma field around the condensate. Our full thermodynamic potential incorporates all cor- 
rections from the medium and vacuum fluctuations in the MS scheme, including logarithmic and 
scale-dependent contributions from quark and sigma bubble-diagrams. Up to 1-loop order, the 
temperature- and chemical-potential-dependent contribution is that of an ideal gas (N c = 3 is the 
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and the quantum vacuum term is given by: 
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In the zero-temperature limit, the medium contribution reduces to: 



ft 



(i) . 

med 






.«.»< -<io g (^fi 



with the effective masses: 



(2.4) 



(2.5) 



(2.6) 



Af, = m q +g(o), 
M 2 = 3(a) 2 -Av 2 



(2.7) 
(2.8) 



The parameters of the lagrangian are chosen such that the effective model reproduces correctly 
the phenomenology of QCD at low energies and in the vacuum, such as the spontaneous (and small 
explicit) breaking of chiral symmetry and experimentally measured meson masses. The conditions 
for fixing the parameters are imposed on the vacuum effective potential, and therefore will be 
modified in the presence of vacuum logarithmic corrections. For more details on the parameter 
fixing procedure, the reader is referred to Ref. [JO]. 

In order to identify the role played by thermal effects and vacuum logarithmic corrections, we 



compare [12] three cases for the effective potential V e ff(a) = U (a) +£2 ren (a): 



(a) including thermal corrections (with T = 30 MeV): £2 ren = Q. } 

(b) considering quantum vacuum terms: Q. ren = Q.^ + ft^g d ; 



(i) . 

med.Th ' 



(c) cold and dense LSMq: £l ren = H 



(i) 

med' 



3. Results 



Let us now analyze the results for the surface tension in the LSMq. For more details, plots, 



discussion and results for other nucleation parameters the reader is referred to Ref. [12]. 

For the zero-temperature chiral transition, the chiral surface tension £ assumes values between 
4 and 13MeV/fm 2 , being well below previous estimates [|l(], 11]. We find that the nucleation pa- 
rameters up to T = lOMeV present variations within ~ 10% of the zero-temperature values. Here, 
motivated by applications to compact stars, we use T = 30 MeV, whenever thermal corrections are 
included. The metastable region in the presence of quantum vacuum terms is ~ 40% bigger than 
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Figure 4: Surface tension for quark-matter nucleation in the cold and dense LSMq (c), including vacuum 
terms (b) and thermal corrections (a). 

in its absence, indicating possible large differences in the dynamics of the system, even though the 
critical chemical potential itself shifts only ~ 2%. 

Fig. H shows a clear competition between vacuum and thermal corrections. As compared to 
the zero-temperature, classical computation (curve (c)), the establishment of an interface between 
quark matter droplets and the hadronic medium costs more with the incorporation of quantum cor- 
rections in the effective potential. On the other hand, finite temperature terms tend to push the 
surface tension down, reaching a minimum of only ~ 3MeV/fm 2 , facilitating the surface forma- 
tion. 

For astrophysical applications, it is sometimes useful to characterize the medium according 
to its density instead of chemical potential. The mapping between density and quark chemical 
potential within the cases we consider is plotted in Fig. g, in which the discontinuities signal the 
first-order phase transition that restores chiral symmetry at high energies. 
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Figure 5: Density in units of the nuclear saturation density (no) as a function of the chemical potential. 



4. Conclusions and perspectives 

We have quantified approximately the homogeneous nucleation in the linear sigma model with 
quarks in the MS scheme at zero and low temperature and finite quark chemical potential, including 
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vacuum and medium fluctuations. All the relevant quantities were computed [jl2p as functions of 
the chemical potential, the key function being the surface tension. 

The inclusion of temperature effects and vacuum logarithmic corrections revealed a clear com- 
petition between these features in characterizing the dynamics of the chiral phase conversion, so 
that if the temperature is low enough the consistent inclusion of vacuum corrections could help 
preventing the nucleation of quark matter during the collapse process. In particular, we show that 
the model predicts a surface tension of ~ 5 — 15MeV/fm 2 , rendering nucleation of quark matter 
possible during the early post-bounce stage of core collapse supernovae. 

The LSMq, however, does not contain essential ingredients to describe nuclear matter, e.g. 
it does not reproduce features such as the saturation density and the binding energy. Therefore, 
the results obtained here should be considered with caution when applied to compact stars or the 
early universe. It is an effective theory for a first-order chiral phase transition in cold and dense 
strongly interacting matter, and allows for a clean calculation of the physical quantities that are 
relevant for homogeneous nucleation in the process of phase conversion. In the spirit of an effective 
model description, our results should be viewed as estimates that indicate that the surface tension 
is reasonably low and falls with baryon density, as one increases the supercompression. First- 
principle calculations in QCD in this domain are probably out of reach in the near future. Therefore, 
estimates within other effective models would be very welcome. 
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